home *** CD-ROM | disk | FTP | other *** search
/ Shareware Overload Trio 2 / Shareware Overload Trio Volume 2 (Chestnut CD-ROM).ISO / dir28 / papc20.zip / REGRESS.PRG < prev    next >
Text File  |  1994-03-24  |  945b  |  45 lines
  1. ; simple linear regression y = a + bx
  2. ; b = (n*sum(xy) - sum(y)*sum(x))/(n*sum(sqr(x)) - sqr((sum(x)))
  3. ; a = mean(y)-b*mean(x)
  4. ; R1  sum(x)
  5. ; R2  sum(y)
  6. ; R3  sum(sqr(x))
  7. ; R4  sum(sqr(y))
  8. ; R5  sum(x*y)
  9. ; R6  n
  10.  
  11. CLRG
  12. LBL dataentry
  13. PROMPT "Enter x and y: "
  14. GTO next  ; only if a null line is entered
  15. 1 6 STO+
  16. DUP 2 STO+
  17. DUP SQR 4 STO+
  18. SWAP
  19. DUP 1 STO+
  20. DUP SQR 3 STO+
  21. * 5 STO+
  22. GTO dataentry
  23.  
  24. LBL next ; continue
  25. ; calculate b = (n*sum(xy) - sum(y)*sum(x))/(n*sum(sqr(x)) - sqr((sum(x)))
  26. 6 RCL 5 RCL * 2 RCL 1 RCL * - 6 RCL 3 RCL * 1 RCL SQR - /
  27. ; calculate a = mean(y)-b*mean(x)
  28. 2 RCL 6 RCL / OVER 1 RCL 6 RCL / * -
  29. PUTS "a = "
  30. 0 PUTX
  31. PUTS "\n"
  32. 8 STO
  33. PUTS "b = "
  34. 0 PUTX
  35. PUTS "\n"
  36. 9 STO
  37. ; calculate RSQ = (a*sum(y) + b*sum(xy) - mean(sqr(sum(y))))/ (sum(sqr(y)) - mean(sqr(sum(y)))
  38. 8 RCL 2 RCL * 9 RCL 5 RCL * + 2 RCL SQR 6 RCL / DUP 0 STO - 4 RCL 0 RCL - /
  39. 4 fix
  40. PUTS "RSQ = "
  41. 0 PUTX
  42. PUTS "\n"
  43. CLST
  44. END
  45.